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Abstract 

The essentially unique reduction of the Euler-Poinsot problem may be per- 
formed in different sets of variables. Action- angle variables are usually pre- 
ferred because of their suitability for approaching perturbed rigid-body mo- 
tion. But they are just one among the variety of sets of canonical coordinates 
that integrate the problem. We present an alternate set of variables that, 
while allowing for similar performances than action-angles in the study of 
perturbed problems, show an important advantage over them: Their trans- 
formation from and to Andoyer variables is given in explicit form. 

Keywords: Euler-Poinsot reduction, Hamilton-Jacobi equation, elliptic 
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1. Introduction 

The Euler-Poinsot problem is a three degrees-of-freedom (DOE) prob- 
lem whose super-int egrable char acter limits the solutions to quasi-periodic 
orbits on two-torus (Fasso, 20051 ). Because of the symmetry with respect to 
rotations about the angular moment um vector , the problem is formulated 
as a 1-DOE Hamiltonian when using lAndoveil variables. Then, the 

essentially unique complete reduction that provides the integration of the 
problem can be performed in different variables, and it is usually done by 
solving the Hamilton-Jacobi equation. 

In the study of the rigid body rotation under external torques the use of 
suitable variables reveals crucial to the solution by perturbat io n methods . 
A com mon trend is to use action- angle variables ( Sadov . 1970a b: Kinoshital . 
1972 ). but other variables can be used instead ( Hitzl and Breakwell 1971 ). 
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The Hamilton-Jacobi equation of the Euler-Poinsot problem in Andoyer 
variabl es can be solved formal ly, without need of specifying the new Hamil- 
tonian ( Ferrer and Lara . 20ld ). We show how this new, unspecified Hamil- 



tonian can be cast into a standard form in which the modulus of the elliptic 
integrals that appear in the solution of the transformation, remains as an 
undetermined state function of the new momenta. Under certain conditions 
imposed to the formal transformation, the modulus is determined by solv- 
ing a system of partial differential equations. The condition we require is 
"simplification" and find a new se t of va r iables that while showing similar 
performances than action-ang les (jSadovl . Il970al lbh. has the benefit of not 



requiring the computation of implicit functions. In addition, we demon- 
strate that Sadov's transformation is also a member of the general family of 
Euler-Poinsot transformations to Andoyer variables. 

2. Complete Reduction of the Euler-Poinsot Problem 

The Hamiltonian of the torque-free rotation is dPepritJ . 



n = [sin^ v/A + cos^ v/B^j (M^ - N')/2 + N'/2C, (1) 

where A, B, and C are the principal moments of inertia of the body, and 
the Andoyer variables are defined by three pairs of conjugate variables: the 
rotation angle on the equatorial plane of the body v and the projection of 
the angular momentum vector on the body axis of maxima inertia A^, the 
precession angle on the invariant plane // and the modulus of the angular 
momentum vector M, and the node angle on the inertial plane A and the 
projection of the angular momentum vector on the axis perpendicular to 
the inertial plane A. Because A, A and fi are cyclic A = Aq, A = Aq, and 
M = Mq are constant, and Eq. ([l]) is a Hamiltonian of 1-DOF. 

The integration may be done by complete reduction. To this goal, we 
look for canonical transformations Tic ■ (A, fi, u, A, Af , N) — )• (£, /i, L, G, H) 
that convert Eq. ([1]) in a new Hamiltonian /C that depends only on momenta. 
Because of the two-torus topology of the Euler-Poinsot problem only two 
momenta are required in /C, and in view of neither A nor A appear in Eq. 
([TJ, we choose h = X, H = A and JC = JC{L, G). 

2.1. Formal Solution of the Hamilton-, J acobi Equation 

In the Hamilton-Jacobi approach, the transformation T/c is derived from 
a generating function in mixed variables S = S{p, v, L, G) such that 

{i,g,M,N) = —— ^ (2) 
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Because // is cyclic in Eq. ([TJ , 5 is chosen in separate variables S = G fj, + 
W{u, L,G). Then, from Eq. ([TJ we form the Hamilton- Jacobi equation 



sin^ V 
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where W may be solved by quadrature. Then, calling /3 = L/G, the trans- 
formation Eqs. ([2]) are 
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VQ 5(i/A) 

sin^ v/A + cos^ z//S - 1/A 



Q 



and 



sin^ i^/A + cos'^ u/B-l/C' 



1/A = 2/C/G2 



(4) 

(5) 

(6) 
(7) 

(8) 
(9) 
(10) 



We only discuss the general case A < B < C. As \/Q must be real for 
all zy, we get 5 < A < C and ^G'^/G < K. < \G'^/B, thus constraining 
the motion to rotations about the axis of maxima inertia. The discussion of 
other cases is left to the reader. 

The transformation equations for i and 5, Eqs. (I1D~([S]), depend on the 
integration of the two quadratures in Eq. ([8]). However, as far as K, de- 
pends only on the momenta G and L, the quadratures in Eq. ([5]) can be 
solved without need of specifying the formal dependence of /C on the new 
momenta, thereby givin g rise to a whole family of canonical transformations 
( Ferrer and Lara . 20ld ). 

The closed form solution of Eq. ([8]) relies on well known changes of var- 
iables. Thus, introducing the parameter / > and the function < m < 1 



/ 



G{B-A) 
{G-B)A' 



m 



{G-A){B-A) 
{G-B){A-A)' 



(11) 
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and the auxiliary variable ip defined as 

Vl + / sin V cos V 

cos ip = — , simp = , 

yl + f sin^ 1/ V 1 + / sin^ i' 

then, the quadratures in Eq. ^ are solved to give 



7 



7 



f + m 
AC 



F(ip\m) — n(— /, Ip I m) 



C- A 



F{ip\m), 



where 
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(l + /)(/ + m)// 



/SA(C7-^) (C-A) 



(12) 

(13) 
(14) 

(15) 



AC{C - B){A-A)' 

F{ip\m) is the elliptic integral of the first kind of modulus m and amplitude 
ip, and Il{—f,'ip I m) is the elliptic integral of the third kind of modulus m, 
amplitude ip, and characteristic — /. It is worth mentioning that for the 
def inition of the elliptic integral of third kind we adhere to the convention 



in ( Bvrd and Friedman! . Il97ll ) 



2.2. The Standard Hamiltonian 

From Eqs. (|10p and (jlip we note that /C is characterized by the identity 

— d- ^ 
2A\ C f + m, 

which can be taken as a definition by assuming that m = m{L, G) in Eq. 
(jl6p . Then, Eqs. (I4])-(I7]) are rewritten 



(16) 



t = : —F['ip\m), 



27 / + m dp 



// + 7 



1 



f + m 



m 



f (3 dm\ 



f + m 2 dp J 



F( 



n(-/,^|m) 



N 
M 



1 — jnsin^ 



(17) 

, (18) 

(19) 
(20) 

Transformations in the literature can be obtained as particular cases of 
the general form Eqs. (ll6D -(pO|). Thus, the new Hamiltonian selected by 
Hitzl and Breakwelll ( 197ll ) is the average of the A ndoyer Ha r niltori ian Eq. 
([1]), which is also the interi nediate Hamiltonian of iKinoshital ()1972l ). while 
a previous proposal of ours ( Ferrer and Lara . 20ld ) transforms the Andoyer 
Hamiltonian to the axisymmetric case. 
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3. New variables 

Searching for simplification in Eqs. (|17p and (|18|) we propose to choose 
1 1 + f dm _ 1 / / (3dm\_ 

2-f f + m d/3 ' f + mV" f + m2 8^) ' ^ ' 

Equations (j2ip can be solved for /? = without need of solving any 

partial differential equation. Furthermore, by squaring /? we can express m 
as a function oi L/G 

m = /[(l + /)GVL2-l]. (22) 
Correspondingly, the new Hamiltonian in new variables is 

whose Hessian never vanishes, and which is formally equal to the uniaxial 
case for a new maximum momentum of inertia, say P, such that 1/P = 
l/A + l/C -l/B. 

Then, the direct transformation is 

I = -F{ij\m) (24) 

9 = + + /) (/ + H// [Hi' I - n(-/, V' I m)] (25) 



L = N \J{1 + f)/{l-msiu^^) (26) 
G = M (27) 

where ip is defined in Eq. (jl2p and ?n is easily computed in Andoyer variables 
from its definition in Eq. ([II]) by noting that A = G^/{2IC) = M^/{2n), 
where Ti is given in Eq. 

The inverse transformation requires using the Jacobi amplitude am to 
invert Eq. 1^ 

ip = -am {e\m), (28) 
where m is computed from Eq. (I22p . Then, from Eq. (I12p we get 



VI + /sn(^,m) . cn(£,m) 

cosi/ = , . sm z> = , (29) 

Vl + /sn^(^, m) Vl + /sn^(A"^) 

where sn, cn, dn, stand for the usual Jacobi elliptic functions. Finally, the 
inverse transformation of Eqs. (j24p - (|27p is completed with 

fi = 5 + (1 + /) (G/L) [£ + n(-/, am {£\m) \ m)] , (30) 
N = Ldn{e,m)/^TTf, (31) 
M = G. (32) 
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4. Transformation to action-angle variables 



Note in Eq. (|'28|) that the variable I is 4Er(m)-periodic, with K{m) the 
complete elliptic integral of the first kind. With a view on perturbations of 
the Euler-Poinsot problem, where elliptic functions would be expanded in 
Fourier series, it could be desired that i be 27r-periodic (an angle). 

The new variable 



vr 



F{^\m), 



2K{m) 

will be obtained by requiring to Eq. (jl7p that 

1 1 + / dm 



IT 



27 / + m d(3' 2K{m) 



(33) 



(34) 



where /3' = L'/G'. Equation (|34p is in separate variables and is integrated 
by quadrature to give, up to an integration constant. 



l3' = -J{l + f) if + m)/f 
vr ^ 



n(-/,m) 



m 



f + m 



K{m) 



(35) 



where !!(—/, m) is the complete elliptic integral of the third kind. Equation 
(j35p defines m as implicit function of V /G' . 

Now, we replace /3' in Eq. (llSh by its value from Eq. (I35p to get 



g' = ^, + J{l + f) if + m)/f 



n(-/,m) 

K{m) 



m 



n(-/,V'|m) 



(36) 



Remarkably, Eqs. (pOj) . ([33|) . (1351) and (1361) re cover the original transfor- 
mation to action- angle variables ( Sadov . Il970al lbl) without need of relying 
on their classical definition^ 

Finally, we note that the inverse transformation from action-angles to 
Andoyer variables requires the computation of m from the implicit function 
Eq. 



^In (jSadovl Il970al lb[) . / = and m = X^, Andoyer variables are {h,ip,(l), L,G,G(;) = 
{X,fi,v,A,M,N), and the action-angles are {f,i',h,I,G,L) = {£' , g' , h' , L' ,G' , H'). Be- 
sides, Sadov's au xiliary angle is ^ = — i/'- Note that there is a typo in the definition of 

in Eq. (4) of (ISadovl ligTOal). whic h should be multiplied by A/C Th i s typ o is easily 
traced in Eq. (2.20) of (|Sadovl ll9 " 



Sadovlll970"bh . but it still remains in (|Sadovl . [l983 : IKozIovI . [2OO0I ') . 
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5. Conclusions 



Action-angle variables are not necessarily the better option for dealing 
with perturbed motion. This fact is very well known for perturbed Kep- 
lerian motion where Delaunay variables are used customarily. The same 
happens to rotational motion where action-angle variables have the incon- 
venience of being related to Andoyer variables through implicit relations. 
But the complete reduction of the Euler-Poinsot problem may be achieved 
in a variety of canonical variables. Indeed, we demonstrate that when solving 
the Hamilton- Jacobi equation of the Euler-Poinsot problem in Andoyer vari- 
ables, the new Hamiltonian can be cast into a standard form as a function of 
the modulus of the elliptic integrals required in the solution, a quantity that 
is consubstantial to the problem. Then, the solution of the Hamilton-Hacobi 
equation can be written formally as a function of the modulus and its partial 
derivatives with respect to the new momenta. Solving these partial deriva- 
tives according to certain criteria provides the desired transformation. In our 
case, we require "simplification" and find a new transformation of variables 
that, while having similar characteristics than action-angles variables, does 
not rely on implicit functions. Besides, we show that the transformation to 
action-angle variables pertains also to this general family. 
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